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The 3-dimensional cube is the only periodic, 
connected cubic graph with perfect state transfer 

Simone Severing 

Department of Physics and Astronomy, University College London, WC1E 6BT London, United Kingdom 

There is perfect state transfer between two vertices of a graph, if a single excitation can travel 
with fidelity one between the corresponding sites of a spin system modeled by the graph. When the 
excitation is back at the initial site, for all sites at the same time, the graph is said to be periodic. 
A graph is cubic if each of its vertices has a neighbourhood of size exactly three. We prove that 
the 3-dimensional cube is the only periodic, connected cubic graph with perfect state transfer. We 
conjecture that this is also the only connected cubic graph with perfect state transfer. 



I. INTRODUCTION 



o . 

' A. State transfer 

' 

Let G = (V, E) be a graph with set of vertices V(G) = {1, 2, n} and set of edges E(G) C V(G) x V(G) - {{i, i} : 
i E V(G)}. The order of G is the number of its vertices. 

Let us consider a system of n spin-1/2 quantum particles with unit XY couplings. The space assigned to the entire 

system is (C 2 )®". Each particle is attached to a vertex of G. The coupling between two particles is nonzero if and 
only if the particles are attached to adjacent vertices. The couplings are then specified by the adjacency matrix of the 
graph. The ij-th entry of the adjacency matrix of G is L4(G)]ij = 1 if {i,j} € E{G) and [^4(G)]y = if {i,j} £ E{G). 
We shall work with the XY model. Let Xj and be the Pauli operators acting on the i-th particle. The Hamiltonian 
£3 governing the dynamics of the spin system can be written as Hxy(G) = 2 _1 [A(G)]ij (XX, +YiYj). Let 

1 i ' {|1) = ex, |2), . . . , \n}} be the standard basis of the space C™. A vector |i) indicates the presence of an excitation at 

vertex i only. With respect to the standard basis, the ij-th entry of the Hamiltonian acting on C" is [HxY(G)]i,j = 
2[A(G)]i j. It follows that the Schrodinger evolution of the excitation is practically induced by a unitary matrix of 
the form Ua(i) = e _4j4 ( G )*, where t g R + . We obtain a probability distribution supported by V{G) by performing 
1 a projective measurement on the state \i)j t ) — UG(t)\ipo). For regular graphs, the "practically" has a much larger 
extension, given that, with constant couplings, any kind of interaction has an Hamiltonian proportional to A(G). 

Given two vertices i,j £ V(G), the fidelity at time t between i and j is the function fG(i,j',t) — \{j\UG(t)\i)\. We 
say that there is perfect state transfer (for short, PST) between the particles i and j at time t if fdhj'^t) = 1 [HI- 
We say that G is periodic, with period t,ii fc(h i;t) = 1 Sometime in the physics literature periodic graphs are 
' said to afford perfect revival (see, e.g., [7[). 
7-H . Even if the topic is not directly discussed here, it is worth noticing that in the XY Z model, the Hamiltonian 
restricted to C™ is proportional to the Laplacian matrix of G (see, e.g., [10]). This fact alone is sufficient to distinguish 
different approaches for the two models, when the graphs considered have generic degree sequences. In this work, we 
■ consider regular graphs only. The result obtained is then also valid for the XYZ model. 

The concept of PST has been introduced in @ and [13| . The recent papers @ and 0, even if not reviews, point 
. . . i out a good number of references embracing the more mathematical aspects around the notion. 



B. Diameter 



A graph H = (W,F) is a subgraph of G if W(H) C V(G) and F{H) C E(G). A subgraph H = (W,F) is an 
induced subgraph of G if H is a subgraph of G and, for every two vertices i,j £ V(H), {i,j} G E(H) if and only if 

{i,j}eE(G). 

The degree of a vertex i is the number of edges incident with i. A path of length I from vertex i to vertex j (if 
there is one) is an induced subgraph with I vertices and I — 1 edges, such that i and j have degree one and all other 
vertices in the path have degree two. A graph is said to be connected if every two vertices are in a path. 
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Let Vij(G) be the set of all paths with end- vertices i and j. The length of a path with end- vertices i and j is 
denoted by The (geodesic) distance between two vertices i and j is defined as d(i,j) = min-p. ^ l(i,j). The 

diameter of a connected graph G is defined as dia(G) = maxij£V(G) d{i,j). Informally, the diameter is the longest of 
the shortest paths. Two vertices i,j G V(G) are said to be antipodal if d(i,j) — dia(G). 

C. Order/distance problems for state transfer 

There is a large and growing literature concerned with the mathematics of state transfer on spin systems. Some 
effort towards a classification may be seen as driven by three recurrent, but essentially unstated problems: 

• General graphs: Given DgN, find the graph G = (V, E) with the smallest possible number of vertices njj — 
\V(G)\ such that d(i,j) = D and fc(i,j',t) = 1 for some t G M + . The set of these graphs is denoted by Gd- 

• Fixed degree graphs: Given D,AeN, find the graph G = (V, E) with the smallest possible number of vertices 
n,D,A = \V(G)\ such that (i) the maximum degree of G is A, (ii) d(i,j) = D and fa(i,j',t) = 1 for some t G K+. 
The set of these graphs is denoted by Gd,a- 

• Regular graphs: Given D, k G N, find the graph G = (V, E) with the smallest possible number of vertices 
n<D,k — \y{G)\ such that (i) G is fc-regular, (ii) d(i,j) = D and fo(i,j;t) = 1 for some t G M+. The set of these 
graphs is denoted by k . A graph is k-regular if all of its vertices have degree k. 

The requirement "smallest possible number" could be replaced with "largest possible number" to state the specular 
versions of the problems. 

D. Examples 

• D = 1: Let P 2 = ({1,2}, {{1,2}}) be the path of length one. Then [Up 2 (i)] 1>2 = -isin(f) and 

max(|[f/p 2 (t)] 1 , 2 |)=/p 2 (l ) 2; 7 r/2) = l. 

Hence, P 2 G Q\. 

• D = 2: Let P 3 = ({1,2,3}, {{1,2}, {2,3}}) be the path of length two. Then [U Pa (i)]i ;3 = -sin(t/y/2) 2 and 

max(|[[/ P3 (t)] 1>3 |) = f Ps (l,S;n/V2) = 1. 

Hence, P 3 G Gi- 
In both cases, PST is between antipodal vertices. 

• D = 3: Let P 4 = ({1,2,3,4}, {{1,2}, {2,3}, {3,4}}) be the path of length three. Let a := t/2. Then [U Pl (t)] M = 

sin (a + a\/5) /10 — i sin (a + a\/5) /2 — i^/5 sin (a — a\/5) /10 — i sin (a — a\/5) /2 and 

max (| [C/p 4 (t)] M |) = /p 4 (i,3;2tt/v^) - sin (n/Vs) w 0.986. 

Hence, P^^Gz- 

E. Cartesian products 

The Cartesian product G = G\ x G2 = (V,S) of two graphs Gi = (Vi,i?i) and G 2 = (V2,S 2 ) has set of vertices 
V(G) = V(Gi) x V(G 2 ) and {{i, j}, {fc, /}} G E (G) if (i) i = k and {j,Z} G S(G 2 ) or ^ j = I and {i,fc} G £(Gi). 
Two facts are important: dia(G) = dia(Gi) + dia(G 2 ); if Gi and G 2 arc fc-regular and Z-regular graphs, respectively, 
then G is (fc + Z)-rcgular. 
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The k- dimensional cube is the graph P 2 xfe (in Fig. [T] P 2 x3 ). We have the following four cases: 



[U p xk (£)]i 2* — 



-sin(i) , i£k = 2l, /odd 
sin (t) k , if fc = 21, I even 
isin(*) fc , if k = 21 + 1, I odd 
—i sin (i) fc , if k = 21 + 1, I even. 



In all these cases, 



xk\\ _ 9 fc „„j j:„/n>xfc 
2 ' vvt - r V- 1 2 / I ~ 

" " < 2 m , where to = D,k 



For P x/£ , we have \V(P£ 

m Gm.m' n m,m 



(\[U p *> (tt/2)] 1j2 »|) = / p x* (1,2*^/2) = 1. 
2 fc and dia(P 2 xfc ) = k- dia(P 2 ) = 21og 2 |F(-P 2 xfe )| = k. Thus, for the fc-regular graphs 



max 

tGR+ 



The fc-dimensional generalization of the 3x3 grid is denoted by P 3 xfc (in Fig. Q] the grid P 3 x3 ). There are four 



types of matrix entries for the graph P 3 xfc : 



[U p xk (t)] 13 k 



sin 



sin (t 



(VV2) 



isin (t/v 7 ^) 
^ -isin (t/V2) 



2k 



2k 



2 k 



21, 



if k = 21, I odd 
if k = 21, I even 
if k = 21 + 1, I odd 
if fc = 21 + 1, / even 



Then, 



max 



(|[t/ p x* (7r/V2)]i l3 *|) =f P ^ (l,3 fc ;7 r/x/2) 



For P 3 xfc , we have V(P 3 



xfc\ 



3* and dia(P 3 x,c ) = fc- dia(P 3 ) = 21og 3 |y(P 3 xfc ) 



2k. 



For P 2 xfe and P 3 xfe PST is between antipodal vertices. The parameters related to PST between two vertices i and 
j in these graphs are given in the following tables jllj]: 



pxfe 


k 


|l/(P 2 X fe )| 


d(i,j) 


pxfe 
^3 


fe 


|l/(P 3 xfc )| 






2 


4 


2 




2 


9 


4 




3 


8 


3 




3 


27 


6 




4 


16 


4 




4 


81 


8 




5 


32 


5 




5 


243 


10 



Notice that P 3 xfe is nonregular for every k. 

A square matrix M of size n consisting of unimodular entries |Afjj | = 1 is called a Hadamard matrix if iJi/t = n/, 
where 7 is the identity matrix and t denotes the Hermitian transpose. In a complex Hadamard matrix, Mij G C 26]. 
The matrix 



UpXk 



(4) V2 



1 -% 
-i 1 



is a complex Hadamard matrix. 



F. Statement of the results 



We prove that the 3-dimensional cube is the only periodic, connected cubic graph with PST. Equivalently, 

Theorem. The 3-dimensional cube, P 2 x3 , is the only periodic, connected cubic (3-regular) graph G with two different 
vertices i and j such that fa(i,j;t) = 1, for some t E R + . 

The statement is verified directly in the next section. Conclusions follow. The proof is based on two known results: 
a periodic, connected regular graph is integral; there are only thirteen connected cubic integral graphs. The proof 
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FIG. 1: (L): The graph W . There is PST between vertices 1 and 8. The spectrum of W is not integral. ( C): The 3-dimensional 
cube. There is PST between vertices 1 and 8. These vertices are antipodal and at distance 3. (R): The graph P 3 x2 . Note that 
the graph has vertices of degree two and three. There is PST between the vertices 1 and 9. These vertices are antipodal and 
at distance four. 



is technically easy, but tedious, because it goes through a seemingly unavoidable case by case analysis. Nonetheless, 
establishing the result is also an excuse for a further step into a systematical exploration of periodic quantum dynamics. 
The proof is interspersed with extra information. The broader aim would be to take a picture of periodic quantum 
dynamics on cubic graphs, even if here we do not state any further general result, beyond a crude analytic compilation 
of matrix entries. 

It is an open problem to prove that P 2 x3 is the only connected cubic graph with PST. 
Conjecture. 3 = {P^ 3 } if -D = 3 and 3 = 0, otherwise. 

Periodicity is not necessary for PST. There are examples of regular graphs that are not periodic but have PST [23| . 



II. PROOF OF THE THEOREM 



A. Integral graphs 

The spectrum of a graph G is the multiset {Ai ^ (G), A^^G), A|.™ r '(G)} of the eigenvalues of A(G). The index 
[mi] in Aj m *'(G) denotes the multiplicity of the eigenvalue Aj(G). For example, the spectrum of the complete graph 
on four vertices, ii^-is {3, — l^}. A graph is said to be integral if its eigenvalues are integers. There is basically one 
survey on this area Q . See also [l5j for the relevant terminology and [l[ for a nontrivial upper bounds on the total 
number of integral graphs with n vertices. There are a few general results that establish a relation between PST and 
integral graphs [l7| . The next statement is directly useful to our purposes: 

PO. A connected regular graph is periodic if and only if it is integral (Corollary 2.3 [17]). 

The converse is not necessarily true. This fact can be observed in the examples discussed below. Moreover, integer 
eigenvalues have a weaker role in nonregular graphs. Let us consider, for instance, a graph W with eight vertices 
{1, 2, 8} and set of edges E(W) = {{1, i}, {j, 8} : 2 < i, j < 7; } (see Fig. [1]). The unitary governing the dynamics 
in W is defined by 



[U w (t)]i 



- sin (V3t) 
—i sin 

cos(V3Y) 2 
(5 + cos (2y/3t)) /6 

— sin (v3i) 



if i = 1 and j = 8 or viz. 

if {i,j}eE(W) 
If i=j = 1,8 
if* = i^l,8 
otherwise 



There is PST between vertices 1 and 8, because [Uw (t/2\/3)]i.8 = — sin(37r/2) 2 = —1. At the same time, the 6x6 
submatrix of Uw (i"/2\/3), excluding the first/last row/column, is 2/3 in the diagonal and —1/3 off-diagonal. The 
spectrum of W is {±2^3, O^}. Thus, W is not an integral graph. 

We shall be interested in a special family of graphs. A cubic graph is a 3-regular graph. All cubic integral graphs 
have been classified and explicitly constructed [UHil]. I n particular, 
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PI. There are exactly thirteen connected cubic integral graphs. 

On the light of the statements PO and PI, a proof of the theorem can be obtained via a case by case analysis. 
As we have seen, and this is an already known fact, the 3-dimensional cube, P 2 x3 , has PST between its antipodal 
vertices. These are vertices at distance three. We shall verify that none of the remaining twelve connected, cubic 
integral graphs affords PST. All graphs considered in this section are periodic. 



B. The complete graph K4, the complete bipartite graph ^3,3, and two connected copies of 1(2,3 

A complete graph on n vertices is a graph K n = ({1, n}, E), where E(K n ) — {{i, j} : 1 < i, j < n}. The ij-entries 
of Uk 4 are given by 



(3 cos (t) + cos (3t) + 4i sin (t) 3 ) /4, if i = j; 
cos (t) sin (t) (—i cos (t) — sin (t)) , if i ^ j. 



Then max teR + ([Cjf 4 (*)]»,») = fu Ki (h i\ 7r /2) = 1, for 1 < i < 4; for every pair of vertices i and j, 
max ieR+ ([U Ki (t)]i,j) = fu Ki (*,i;7i"/ 4 ) = V 2 - Note that 



[^(T/4)]i,i 



J_l (l + i)/2, ifi = j; 
x/2 1 -(l + i)/2, if i ^ j. 



gives a complex Hadamard matrix. 

A graph G = (V = Vi U V2, -E) is bipartite if each vertex in Vi is adjacent to vertices in V 2 only and viz. A complete 
bipartite graph is a bipartite graph K q _ p = (V = V\ U V2,E) such that |Vi| = p, | V2 1 = q and _E = {{i, j} : i £ V\ and 
j e V2}. The graph ^3,3 is on six vertices; ^(#3,3) = with \A\ = \B\ = 3. By definition, {i,j} <= E(K^ 3 ) if 

and only if « e A and j e £>. The spectrum of if^ is {±3,0^} and dia^s^) = 2. The ij-entries of Uk 3}3 are as 
follows: 



(2 + cos (3t)) /3, i = j; 

(-l + cos(3t))/3, {i,j}$E(K 3i3 ); 
-isin(3i)/3, {i,j}6£(^). 



Hence, max ieR + (| [t/ir 33 (<)]»,» |) = fu K33 (i,i;2n/3) = 1, for 1 < i < 6. For the off-diagonal entries, we need 
to distinguish two cases: (i) max teK + (| [U K3t3 (t)]i,j |) = /c/ K3 3 (i,j;n/3) = 2/3, if i,j e A or i,j e 5; ^ 
max feI+ ([|*7fc 3 , 3 (*)]»,.?' |) = fu K33 («,i;T/6) = 1/3, if i G A and j e B. When f = tt/2, 



[^s>/2)k; 



2/3, i = i; 

-1/3, {*, j} ^3,3); 
i/3, 



Let DK 2 ,3 be the graph on ten vertices obtained from two disjoint copies of if 2,3, say K 23 and iff, 3; by adding 
three edges between the vertices of degree two in if 23 and iff 3. The spectrum of DK 23 is {±3, ±2, il^, O^ 2 '} and 
di&(DK 2 .3) = 3. The structure of Udk 3 3 consists of various kind of entries: 



UDK 2 ,'A t ) 



Ol 


02 


03 


03 


03 


0,4 


a 4 


a 4 


«5 


«5 


«2 


fli 


03 


«3 


«3 


04 


04 


a4 


a 5 


a 5 


az 


A3 


a 6 


a 7 


a 7 


a 8 


a 9 


a 9 


04 


04 


0,3 


«3 


a 7 


a 6 


a 7 


a 9 


a 8 


a 9 


04 


04 


03 


03 


a 7 


a 7 


a 6 


a 9 


a 9 


a 8 


04 


04 


04 


04 


«8 


09 


a 9 




ai- 


a 7 


«3 


03 


O4 


04 


a 9 


o 8 


a 9 


a 7 


de 


a 7 


«3 


03 


O4 


04 


o 9 


a 9 


a 8 


a 7 a 7 a 6 


«3 


03 


«5 


a 5 


04 


ft4 


04 


03 


03 


03 


Oi 


02 


«5 


a 5 


04 


ft4 


04 


03 


«:i 


03 


fl2 


Ol 
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where 

01 = (5 + 3 cos (2i) + 2 cos (3t)) /10, 

a 2 = (-5 + 3 cos (2t) + 2 cos (3*)) /10, 

a 3 = -i (sin (2f) + sin (3i)) /5, 

a 4 = (— cos (2i) + cos (3t)) /5, 

a 5 = i (3 sin (2t) - 2 sin (3t)) /10, 

o 6 = (10 cos (t) + 2 cos (2t) + 3 cos (3t)) /15, 

a 7 = (-5 cos (t) + 2 cos (2t) + 3 cos (3f)) /15, 

a 8 = -« (10 sin (t) - 2 sin (2f) + 3 sin (3f)) /15, 

a 9 = i (10 sin (i) + 2 sin (2b) - 3 sin (3i)) /15. 

By considering ai and etg, we can see that max tgR + (| [Udk 2 3 (t)]i,i |) — Jdk 23 (*j*;27t) = 1, for every i. However, 
maxi<j<9 ; jjti t e max tel + ( | a j | ) = (5 — 5 (—1 — v5) /4) /10 ~ 0.9. The maximum is attained by a 2 for t = 2n/5. The 
segments in the matrix U dk 3 3 (t) help visualizing its structure and these do not have a mathematical meaning. We 
shall make a consistent use of this graphic tool also in the next cases. The graphs K3 3 and DKi$ are in Fig. ([2]). 



C. The graphs C 3 + K 2 and C 6 + K 2 



The (Cartesian) sum G = G\ + G2 = (V,E) of two graphs G\ — (V\,Ei) and G2 = (^2,^2) has set of vertices 
V(G) = V{Gi) x V(G 2 ) and {{i, j}, {k, I}} £ E (G) if (i) {i,k} G E(G X ) or fiij {j,?} e ^(G 2 ) [12. We denote by 
C„ the n-cycle: V(C n ) — {1, 2, n} and {i, (i + 1) modn} € E(C n ). 

The graph C3 + has two cycles of length three and it can be drawn as a prism with triangular basis. It is 
the undirected version the Cayley digraph of the dihedral group Dq generated by the standard set. Its spectrum is 
{3, 1, Ol 2 !, — 2 1 2 !}. Given the symmetry, the unitary matrix has a neat structure: 



Uc 3 +K 2 (t) — 



ai 


«2 


«2 


03 


04 


«4 


02 


ai 


a 2 


04 


03 


04 


a 2 




01 


04 


04 


O3 


03 


«4 


04 


ai 


02 


O2 


04 


03 


04 


02 


Ol 


a 2 


04 


O4 


03 


02 


«2 


01 



with 



-it 



2e 



2it 



ai = [2 + e 
a 2 = -ie~ lt / 2 (sin(i/2) 
a 3 = (2 - e- lt - 2e 2it 4 
a 4 = (-1 - e- 1 * + e 2rf - 



e- 3it ) /6, 
+ sin(5t/2))/3, 
e" 3lt ) /6, 

,-3it\ 



e 



/6. 



From this, max 2 <j<4 max t6R + (\aj\) w 0.9, for j ; = 3 and t = v3 + 7r/17. 

The graph G@ + A" 2 is on twelve vertices. It has two cycles of length six and it can be drawn as a prism with 
hexagonal basis. In fact, in analogy with C 3 + K2, it is the undirected version the Cayley digraph of the dihedral 
group Di 2 generated by the standard set. As we have done for the other cases, we explicitly write down the unitary 
matrix. Let 



.4 



«1 


«2 


03 


O4 


03 


«2 


(12 


a 1 


2 


O3 


04 


03 


03 


«2 


Oi 


a 2 


03 


04 


O4 


a 3 


02 


Ol 


02 


03 


03 


04 


O3 


02 


Ol 


02 


02 


03 


O4 


03 


2 


Ol 



and 



05 


«6 


07 


08 


07 


06 


06 


a 5 


06 


7 


a 8 


07 


a 7 


a 6 


a 5 


a 6 


07 


a 8 


o 8 


a 7 


06 


a 5 


a 6 


a 7 


a 7 


as 


07 


(7,5 


05 


06 


06 


a 7 


Og 


07 


06 


a 5 



Then 



Uc e +K 2 (t) 



A B 
B A 
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FIG. 2: From the left: The complete bipartite graph ^3,3. The graph obtained by connecting together two copies of ^2,3. The 
graph C3 + Ki. The graph Ca + Ki. Although periodic, none of these graphs has PST. 



and 



For j ^ 1, 



ai = (2 + cos (t) + 2 cos (2t) + cos (3f)) /6, 
a 2 = -i (sin (i) + sin (2t) + sin (St)) /6, 
a 3 = (-1 + cos (i) - cos (2t) + cos (3*)) /6, 
a 4 = i (sin (i) - 2 sin (2t) + sin (3i)) /6, 
0.5 = 1 (sin (t) — 2 sin (2t) — sin (3t)) /6, 
a 6 = - (1 + 2 cos(t)) sin (i) 2 /3, 
07 = — i (sin (i) + sin (2t) — sin (3t)) /6, 
a 8 = 16 (cos(t/2) 2 sin(t/2) 4 ) /3. 



max (I a,- 1) < 

tm+ J 



21/50, j = 2 and t w tt/6; 

= 2/3, j '■ — 3 and t = 7r; 

29/50, j = 4 and f = 77/20 
« 9/20, j = 5 and t = 19/10 

« 1/2, j = 6 and t w 1.21 
« 9/25, j = 7 and i 1.35 

64/81, j = 8 and t « 1.9 



The graphs 63 + ^2, and 6*6 + if 2 are represented in Fig. @. 



D. The Petersen graph, a graph on ten vertices, and L (S (K4)) 

The Petersen graph, P, illustrated in Fig. (|3J), is one of the best-studied single objects in the graph-theoretic 
literature. The Petersen graph has two cycles of length five. It is vertex-transitive but not a Cayley graph, with 
spectrum {3, 1",— 2"} and dia(P) = 3. It is strongly regular with parameters (10,3,0,1). The symmetry appears 
also in Up(t), which reflects faithfully the structure of P: 



U P (t) = 



a 1 


«2 


«3 


0-3 


a 2 


a 2 


a3 


a3 


a3 


a3 


a 2 


Oi 


«2 


a3 


a3 


a3 


a-2 


a3 


a3 


a3 


03 


2 


ai 


a 2 


a3 


a3 


a3 


«2 


a3 


a3 


A3 


a 3 


a-2 


ai 


a 2 




a3 


03 


(12 


0.3 


ai 


03 


a 3 


a 2 


a\ 


a 3 


03 


0.3 


03 


a 2 


a 2 


a 3 


a 3 


«3 


a3 


a x 


a3 


a 2 


a 2 


03 


0*3 


«2 


«;s 


a 3 


a3 


a3 


a 1 


a3 


a-2 


«2 


G3 


a 3 


a-2 


a 3 


a3 


a 2 


a3 


ai 


a3 


«2 


G3 


a 3 


a$ 


«2 


a3 


a 2 


«2 


a3 


ai 


03 


A3 


«3 


as 


a3 


a 2 


a3 


«2 


(12 


«3 


a\ 
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where 

d! = e- 3lt (1 + 5e 2lt + 4e 5lt ) /10, 
a 2 = e- 3it (3 + 5e 2lt - 8e 5l *j /30, 
a 3 = e- 3it (3 - 5e 2lt + 2e 5lt ) /30. 

Thus, max :)= 2 ! 3 max teR + (\cij\) = 8/15, for j = 2 and t = 7r. More generally, [£/p(7r)]i.j = —1/5 if z = j; 2/15 if 
{h j} i E{P) and -8/15, otherwise. 

There is another cubic integral graph on ten vertices, obtained by replacing with triangles two nonadjacent vertices 
ofX 3 , 3 . Denoted by Z, it has spectrum {3,2, ll 3 l,-l -2 1 3 !}. Fig. ^ contains a drawing. The unitary matrix 
obtained from Z is 



Uz(t) = 
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a 5 


a 7 


04 


as 


«2 


as 


as 


a 7 


a 6 
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a 6 
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04 


a 6 


a 5 



The entries are 
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= e 


-Sit 


(1 + 5e 2it + 4e 5it ) /10, 
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-Sit 


(3 + 5e 2it - 8e 5it ) /30, 
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= e 


-Sit 


(3 - 5e 2it + 2 5it ) /30, 


a 4 


= e 


-Sit 


{e il - if (3 + e il + 4e 2it + 7e 3it ) /30, 


a 5 


= e 
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(3 + 5e lt + 5e 2it + 10e 4it + 7e 5it ) /30, 


06 


= e 


-Sit 


(3 + 5e i4 - 5e 4it - 3e 5it ) /30, 


«7 


= e 


-Sit 


(3 - 5e lt + 5e 4it - 3 5lt ) /30. 



We then obtain maxixjxgy^i^ max teK + (\aj\) ~ 0.85, for j = 4 and twi - 5/6. Note that 01, a 2 , and 03 give the 
same dynamics as the Petersen graph. 

The graph X = L (S (K4)) is constructed by replacing each vertex of K4 with a triangle (see Fig. ©). The triangles 
are then connected by independent edges. The notation indicates the line graph of the K4 subdivision. Its spectrum 
is {3, ±2t 3 l, O' 2 ! , — ll 3 '} and dia(X) = 3. The unitary has some symmetry: 
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FIG. 3: (L): The Petersen graph. (C): The graph Z on ten vertices. (R): The graph L(S(Ki)) on twelve vertices. These 
graphs are periodic without PST. 
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For j ^ 1, 



01 = (2 



0.2 
«3 



V 3e lt + 3e- 2it 
-3« (_2 - 5e i4 



3e 2 



/12, 



2e 3lt + 2e 4lt + 3e 5lt ) /24, 



-3ii 



1 



it 



2e 



3it 



_4it 



a 4 = e" 3lt (e i4 



(e i4 - l) (2 + 3e i4 + 4e 2it + 3e 3it ) /24 



nit 



)/12, 



a 5 
a 6 



-3<i 



max (|a, |) < 



-1) (2 4 
l) 3 (e lt 



«l/2, 
«l/2, 
«l/4, 
= 2/3, 



3e lt + 3e 2it ) /24, 
+ 1)/12. 



j = 2 and f « 7r/4 
j = 3 and £ « 7r/4 
j = 4 and t ~ 6/5 
j = 5 and £ = 7r 



V3/4, j = 6 and i = 2tt/3. 



E. The Desargues graph and its cospectral mate 

The bipartite double cover of the Petersen graph is called Desargues graph. There are many different notations for 
this graph. We adopt H 5t2 - The Desargues graph is on twenty vertices and its spectrum is {±3,±2M,±]J 5 ]}. The 
graph i?5,2 has a cospectral mate, which we will denote by H' 5 2 - This is a nonisomorphic graph with the same spectrum. 
Two graphs G and H are said to isomorphic if there is a permutation matrix Q such that QA(G)Q T = A(H). It is 
clear that two isomorphic graphs have the same spectrum. The converse is not necessarily true. Indeed, H^^ and 
H' 5 2 are a counterexample. The graphs H§ 2 and H' b 2 are drawn in Fig. (0}. Let us define the arrays 
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Let M be the matrix obtained by permuting the lines (rows and columns) of a square matrix M such that line i in 
M is line n — j + 1 in M. One can observe that 



where 



A B 
B A 



a 1 
02 
03 

r/,4 
05 
rtfi 
07 



e -3it ^ 

e" 34 * (3- 



4e 4 * 
2e 4 * 

-I) 4 



4e 



5it 



/20, 



_ „-3it 



3 + 8f 



- 5e 2lt + 5e 4lt 

- 5e 24 * - 5e 4lt + 2e 5lt + 3e 6lt ) /60 
(3 + 4e 4 * + 3e 2it ) /60, 

- 5e 2tt - 5e 4lt - 8e 5it 



_ e -3it 



— e 

—e 



Sit ^it 
-3it ^it 



-Sit 



Ae 2U 

l) 3 (3 + 2e lt + 2e 2lt 
5 (l + e")/20, 

7e lt + 7e 2lt 



1 1 + 4e' 



1 



3e 6it ) /60, 

3ifj 



<■-') /20, 



D 3 (3 



/60, 
3e 3lt ) /60. 



From these functions, we can see that maxi<j<g ; j^i max t6R + (\aj\) s» 0.83, for j = 7 and t s» 575/250. When i = 7r, 
the probability amplitude is supported by all vertices in a class of the bipartition; in other words, the matrix Ujj' 2 ("") 
is block-diagonal with two 10 x 10 blocks corresponding to the classes. The functions in the above equations completely 
specify the dynamics in H^^, since H' 5 2 and H' 5 2 are cospectral. One can verify that the matrix Uh 5 2 (t) is obtained 
by rearranging the entries of Ujj' (t) . 



F. The Nauru graph and the Tutte-Coxeter graph 



The Nauru graph, A/24, is the only cubic symmetric (arc- transitive) graph on 24 vertices. It is bipartite, with 
spectrum {±3, ±2l 6 l, il^ , 0^ } and dia(JV24) = 4. A brief parenthesis: the Foster census of cubic symmetric graphs 
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highlights that a large portion of periodic cubic graph is symmetric [141 ] . Clearly, the two sets do not coincide. There 
are exactly seven different kind of entries in Ujy 2i (t). This can be seen as a 2 x 2 block matrix. The blocks (1,1) and 
(2, 1) are below. The other blocks are just their rearrangements: 
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where 

ai = (2 + 3 cos (t) + 6 cos (2t) + cos (3t)) /12, 
a 2 = - ((l + 2cos(i))sin(i) 2 ) /6, 

a 3 = 8 (cos (t/2) 2 sin (i/2) 4 ) /3, 

a 4 = 2 ((1 + 2 cos (i)) sin (t/2) 4 ) /3, 

a 5 = -i (sin (i) + 4 sin (2i) + sin (3t)) /12, 
a,6 = i sin (i) 3 /3, 

a 7 = -i (sin (t) - 2 sin (2i) + sin (3i)) /12 
are all the the different entries of the unitary matrix. From these functions, we can write 



max (|Oj|) = < 



= 1, 


j = 1 and t = 27T 


«l/2, 


j = 2 and i w 7r/4 


«l/4, 


j = 3 and t s» 6/5 


< = 2/3, 


j = 4 and t = tt; 


« 9/20, 


j = 5 and t « 7r/4 


= 1A 


j = 6 and t = tt/2 


k «3/10, 


j = 7 and i w 3rr/4 



The Tutte-Coxeter graph, TC, is illustrated in Fig. ([5J. This is the largest cubic graph giving a periodic dynamics. 
Its eigenvalues are {±3, ±2^1 , QP°J}. Let us describe how to write down the unitary matrix Utc- In Fig- © we 
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have ordered the vertices anticlockwise. Each even vertex is connected to odd vertices, and viz. As we have seen in 
previous examples, TC is bipartite and its unitary has a block structure: 



Utc 



A B 
B C 



The ij-th entry of the block B is 



-i(6sin(2i) + sin(3t)) /15, {i,j} e £(TC); 
i (3 sin (2i) - 2 sin (3t)) /30, {i, j} </ £(TC). 



The index i runs over the odd numbers 1,3, ...,29; j over the even ones, 2,4, ...,30. For any chosen t, none of the 
entries of B can have unit absolute value. The blocks A and C have the same entries but rearranged: 
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where 



ai = (5 + 9 cos (2t) + cos (3i)) /15, 
a 2 = (-5 + 3 cos (2t) + cos (3i)) /30, 
a 3 = 2 ((7 + 8 cos (i)) sin (i/2) 4 ) /15, 



and 



max ( 

teR+ 



= 1, j = 1 and t = 2tt; 
s3/10, j = 2andtwl/4; 
13/50, j = 3 and f w 91/50. 



This last equation concludes the proof of the theorem. Notice that Utc(k) is 2 x 2 block diagonal: the diagonal 
entries are —13/15; the off-diagonal ones 2/15. The two blocks are Grover matrices. When t = ir/2, all entries in the 
off-diagonal blocks are equal to i/15. 



III. CONCLUSIONS 



We have shown that the 3-dimensional cube is the only periodic, connected cubic graph with PST. The proof is 
based on known results of graph theory and basic definitions of quantum dynamics on graphs. Our proof goes through 
the relevant cases, which can be interpreted as a systematical exploration of periodic quantum dynamics on cubic 
graphs. A necessary and sufficient condition for PST would give a shorter, more elegant proof. However, any method 
for the same task requires the different cases. There is a small plus in the approach adopted here: we have written 
down explicitly the unitary matrices that specify the dynamics. The matrices are potentially useful in further work. 
The well-known mathematical scenario of state transfer together with various facts observed in this paper suggest 
some reflections. These may be worth a mention. 
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FIG. 5: (L): The Nauru graph. (R): The Tutte-Coxeter graph. This is the largest periodic cubic graph. The Nauru graph and 
the Tutte-Coxeter graph do not have PST. 



A. Extremality 

The literature on interconnection networks for telecommunications, parallel computers and distributed systems 
contain many optimization scenarios dealing with order/degree/diameter (l6| . The most famous one is perhaps the 
degree/diameter problem [20| : given natural numbers D and A, find the largest possible number of vertices 7t.a,d in a 
graph of maximum degree A and diameter < D. The graphs achieving the upper bound are known as Moore graphs. 
Interestingly, in the setting of PST, the original problem was somehow the opposite one [lj|. In fact, the main point 
does not seek to maximize the number of vertices without compromising reachability, but to perform long distance 
communication, with a minimum of physical resources. The number of particles equals the order of the graph and it 
is therefore a resource. The order/distance problems for state transfer are then a new scenario, where graphs with 
extremal properties can be defined via some parameter associated to a dynamical process, instead of a topological 
condition. The landscape becomes even richer, if we consider more than a single excitation, and lift the analysis to 
the graph powers defined in [4|. These considerations remind that the problem of identifying extremal graphs with 
respect to state transfer is still a mostly unexplored area of research. 

B. A notion of persistency 

As a generalization of the rule inducing a continuous-time random walk, the operator Ua(t) = e~ %A( fi) t has been 
subject of intense study (see pTj and the references therein). In analogy with random walks, given the importance 
of these objects for constructing distributions, a principal direction has pointed uniform sampling. Additionally, still 
from the quantitative side, it is useful to deal with problems related to single vertices, or pairs of vertices, more than 
the entire graph. For example, an instance of relevant parameters is a quantum version of the hitting time. We would 
like to propose a notion that intends to quantify how the fidelity for state transfer between two vertices is constant 
in an interval, including small fluctuations. In the previous section, we have seen that half of the diagonal entries 
of U C6 +K 2 (t) have the form ax = (2 + cos (t) + 2 cos (2i) + cos (3t)) /6. When t £ {tt/2, 3tt/2}, a x = 0. Roughly 
between 7r/2 and 37r/2 there is a plateau: the value of the function is about 1/3; it is exactly 1/3, for t = tt. Starting 
the process from some vertex i, if at time t = tt we sample from the distribution, we are going to obtain i with 
probability 1/9. This probability is not spicked, but fairly stable around tt. Given a graph G, the e- persistency of 
a pair of vertices i,j £ V(G) is the length of the longest interval T C [0,27r] such that there exists a value k for 
which k — e < \[Ua(t)]ij\ < k + e, with e > 0, for every t € T. By looking at all vertices, one could define the 
maximum persistency or the average persistency, if extending the definition in the obvious ways. A process with 
higher persistency requires less clock precision for sampling. A first sight reason for giving attention to persistency 
could arise from a connection with energy transfer problems (l2|. Let us add an additional remark. Decoherence 
has been shown to behave as a natural smoothing mechanism on probability amplitudes (l9j . Because of this fact, 
does decoherence increase persistency? With a similar scope but a on a different line, is the presence of decoherence 
compatible with PST at all? 
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C. Discrete probability transfer 



The evolution governed by the Hamiltonian Hxy(G) is driven by the exponentially smaller licit), when taking 
into account a single excitation only. While A{G) faithfully represents G, the matrix licit) does not preserve its 
topological structure. In other words, [Ua(t)]ij 7^ does not imply [A(G)] ii:) = 1. A discrete evolution on G could be 
defined (in some cases 525]) by a unitary Wg, with [Wig]t,j ^ if and only if [A(G)]ij — 1. Such a process does not 
describe the transfer of a single excitation, but it only allows to create a probability distribution supported by V(G). 
The process is discrete and it follows the iteration WcW^T 1 i — > W G . If there is t £ N such that |(j|Wg|i)| = 1 then 
we have perfect probability transfer from vertex i to vertex j. Unless we make use of some kind of lifting [2J (e.g., the 
introduction of extra degrees of freedom), this is the closest analogue to the continuous object licit), even if this one 
does not always exist. For example, we can construct on K4 the unitary below: 

"0-11 1 



1110 

However, it is simple to observe that each power of Wka has one of the two zero-patterns 
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where * denotes a generic nonzero entry. This is sufficient to show that Wk a does not give perfect probability transfer 
in K4. In more complicated situations, the zero-patterns of matrix powers are not immediately available to imply 
a general statement. To verify that a graph G does not enjoy the property, we should study the spectra of unitary 
matrices with the same zero-pattern of Wc- Since the problem involves both spectra, zero-pattern, and an optimization 
procedure, its flavour reminds of the matrix analysis questions approached in Q or various parametrizations coming 
from graph matrices |18j . In our context, the use of semidefinite programming techniques does not seem immediately 
useful, because the matrices are not stochastic, but in fact unitary. 
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